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We refine the relativistic constituent quark model developed in our previous papers to include 
' the confinement of quarks. It is done, first, by introducing the scale integration in the space of a- 

\ parameters, and, second, by cutting this scale integration on the upper limit which corresponds to 

an infrared cutoff. In this manner one removes all possible thresholds presented in the initial quark 
' diagram. The cutoff parameter is taken to be the same for all physical processes. We adjust other 

model parameters by fitting the calculated quantities of the basic physical processes to available 
experimental data. As an application, we calculate the electromagnetic form factors of the pion and 
the transition form factors of the cj and r) Dalitz decays. 
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I. INTRODUCTION 

m 

> . 

' The analysis of the structure of hadronic matter is one of the key tasks of modern particle physics. The generally 
T-H I accepted view is that hadrons are made up of the quark and gluon degrees of freedom of Quantum Chromo Dynamics 
■ (QCD) (see, e.g. the monograph An attempt to give a quantitative field theoretic description of hadrons and 
I their interactions with quarks and gluons has to go beyond perturbation theory and necessarily implicates the use 
^NJ ' of nonperturbative methods. Such phenomena as hadronization, i.e. how hadrons are constructed from quarks, and 
confinement, the empirical fact that quarks have not been detected in isolation, can only be understood via non- 
0^ . perturbative methods. Significant progress has been made in constructing various quark models of hadrons which 
implement different features of QCD. For example, potential models provide simple tools that allow one to describe 
the hadron spectrum. However, the use of quantum mechanical potential models cannot explain the confinement 
of light quarks because quark creation and annihilation effects are essentially nonperturbative. In the approaches 
^ , based on quantum field theory one can understand quark confinement as the absence of quark poles and thresholds 
' in Green's functions and matrix elements. 

In the Quark Confinement Model (QCM) Ref. [2] quark confinement was implemented by assuming that, at low 
energies, the constituent quark interacts with some given vacuum gluon configurations. As a result the quark has no 
fixed pole mass and there are no quark poles in the Green's functions and matrix elements. Models based on results 
obtained via QCD's Dyson-Schwinger equations (DSEs) 3] possess the feature that quark propagation is described by 
fully dressed Schwinger functions. Dressing of the quarks eliminates the threshold problem and thus one is effectively 
working with confined quarks. Using the framework of quantum field theory a promising approach for the description 
of composite particles as bound states of their constituents was suggested in [4]. 

Over the past years we have developed a relativistic constituent quark model and have applied the model to a 
large number of elementary particle processes Q-Q- The relativistic constituent quark model can be viewed as an 
effective quantum field approach to hadronic interactions based on an interaction Lagrangian of hadrons interacting 
with their constituent quarks. The coupling strength of the hadrons with the constituent quarks is determined by 
the compositeness condition Zh = Q where Zjf is the wave function renormalization constant of the hadron. The 
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hadron field renormalization constant Zn characterizes the overlap between the bare hadron field and the bound 
state formed from the constituents. Once this constant is set to zero, the dynamics of hadron interactions is fully 
described by constituent quarks in quark loop diagrams with local constituent quark propagators. Matrix elements 
are generated by a set of quark loop diagrams according to an expansion. The ultraviolet divergences of the 

quark loops are regularized by including vertex form factors for the hadron-quark vertices which, in addition, describe 
finite size effects due to the non-pointlike structure of hadrons. The relativistic constituent quark model contains only 
a few model parameters: the light and heavy constituent quark masses and the size parameters that describe the size 
of the distribution of the constituent quarks inside the hadron. 

In the light quark sector this approach was successfully applied to describe the electromagnetic properties of the 
pion and the sigma meson the electromagnetic form factors and magnetic moments of nucleons, and semileptonic 
decays of the light ground state baryon octet In the heavy quark sector we have calculated the baryonic Isgur- 
Wise functions, decay rates and asymmetry parameters for semileptonic decays of baryons, non-leptonic decay rates, 
and one-pion and one-photon transitions of heavy fiavored baryons M] . This technique was also applied to study the 
semileptonic decays of the double heavy baryons and the Sc-meson Q ■ The same bound-state formalism was recently 
applied to the analysis of exotic hadron states -so-called hadronic molecules- in which the constituents are hadrons 
themselves rather than quarks as in the present approach [9]. 

The local form of the quark propagators used in the relativistic constituent quark model can lead to the appearance 
of threshold singularities corresponding to free quark production in transition amplitudes. As a result, applications of 
the relativistic quark model had to be restricted to ground state mesons and baryons with masses less than the sum of 
the constituent quark masses, and processes with relatively small energies. This poses problems for the description of 
e.g. light vector mesons, e.g. (p, if*), and excited states where the particle mass exceeds the sum of the constituent 
quarks. 

In the present work we propose a refinement of our previous quark model approach by effectively implementing 
quark confinement into the model. In the relativistic constituent quark model matrix elements are represented by a 
set of quark loop diagrams which are described by a convolution of the local quark propagators and vertex functions. 
By using Schwinger's a-representation for each local quark propagator and integrating out the loop momenta, one 
can write the resulting matrix element expression as an integral which includes integrations over a simplex of the 
a-parameters and an integration over a scale variable extending from zero to infinity. By introducing an infrared 
cutoff on the upper limit of the scale integration one can avoid the appearance of singularities in any matrix element. 
The new infrared cutoff parameter A will be taken to have a common value for all processes. We determine the 
parameters of the model by a fit to available experimental data. As a first step, we apply our approach to evaluate 
the electromagnetic form factors of the pion and the transition form factors of the Dalitz decays P — )■ 7 and 
V Pl+l'. 

The paper is structured as follows. In Sec. |TT1 we discuss the basic notions of the relativistic constituent quark 
model including a discussion on how gauge invariance is implemented in the model. We then go on to explain how 
to effectively implement infrared confinement. In Sec. IIIII we determine the model parameters by a fit to the decay 
constants of the tt and p mesons and then apply the covariant constituent quark model to evaluate the decays tt" — 7 7 
and p ^ TTj including also the form factor behavior of the electromagnetic transitions. In Sec. IIVI we extend our 
approach to the strange, charm and bottom sectors. We calculate the leptonic decay constants of both pseudoscalar 
and vector mesons, and the electromagnetic and leptonic decay widths. We calculate transition form factors and the 
widths of their Dalitz decays, and compare the results with recent experimental data. Finally, in Sec.|Vl we summarize 
our results. In Appendix A we describe in some technical detail how the one- loop integrations involving tensor loop 
momenta have been done. In Appendix B we proof by explicit calculation that the p — >■ 7 transition matrix element 
discussed in the main text is gauge invariant. 



II. THEORETICAL FRAMEWORK 



A. Lagrangian 

The relativistic constituent quark model is based on an effective interaction Lagrangian describing the coupling of 
hadrons to their constituent quarks. In this paper we limit ourselves to the meson sector in which mesons are described 
as quark-antiquark bound states. An extension of the model to baryons (three-quark states) and multiquark states 
is straightforward. The coupling of a meson M[qiq2) to its constituent quarks qi and 52 is described by the nonlocal 
Lagrangian 

^t^t{x) = 9mM{x) dxi dx2FM{x,xi,X2)qi{xi)TM>^Mq2{x2) + h.c. (1) 
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Here, Xm and Vm are Gell-Mann and Dirac matrices (or a string of Dirac matrices) chosen appropriately to describe 
the flavor and spin quantum numbers of the meson field M{x). In the case of pseudoscalar mesons we introduce 
singlet-octet mixing with a mixing angle of Op = —18°, while the vector mesons are assumed to be ideally mixed. 
The vertex function Fm^x, xi,X2) characterizes the finite size of the meson. To satisfy translational invariance the 
vertex function has to obey the identity Fm{x + a,xi + a^X2 + a) = Fm{x, xi,X2) for any given four-vector a. In 
the following we use a specific form for the vertex function which satisfies the above translation invariance relation. 
One has 

Fm{x, Xi,X2) = 5^^\x-^W^Xi)^M{{Xl-X2f\ (2) 

where $m is the correlation function of the two constituent quarks with masses mi and TO2- The variable Wi is defined 
by Wi = mi /(mi + 7712) so that wi -I- u'2 = 1. In principle, the Fourier transform of the correlation function, which 
we denote by $m(— ^^), can be calculated from the solutions of the Bethe-Salpeter equation for the meson bound 
states 3. In Refs. [5] it was found that, using various forms for the vertex function, the basic hadron observables are 
insensitive to the details of the functional form of the hadron-quark vertex form factor. We will use this observation 
as a guiding principle and choose a simple Gaussian form for the vertex function <I>m(^^^)- The minus sign in the 
argument of <I'a/(— is chosen to emphasize that we are working in Minkowski space. One has 

^M{-l')^exp{iyAlj) (3) 

where the parameter Am characterizes the size of the meson. Since P turns into —P in Euclidean space the form Q 
has the appropriate fall-off behavior in the Euclidean region. We stress again that any choice for <I> m is appropriate 
as long as it falls off sufficiently fast in the ultraviolet region of Euclidean space to render the Feynman diagrams 
ultraviolet finite. 

In the evaluation of the quark-loop diagrams we use the free local fermion propagator for the constituent quark 

S,ik) = (4) 

mq— f — l€ 

with an effective constituent quark mass mq . The local form of the quark propagator will lead to the appearance of 
threshold singularities corresponding to free quark production. It is for this reason that we had restricted previous 
applications of the relativistic quark model to the lowest- lying states that satisfy the condition m^^ +TOg2 > ™m, and 
to processes with relatively low energies. 

The coupling constant gM in Eq. ([1]) is determined by the so-called compositeness condition. The compositeness 
condition requires that the renormalization constant Zm of the elementary meson field M{x) is set to zero, i.e. 



where (mj^,j) is the derivative of the mass operator corresponding to the self-energy diagram Fig. [T] One has 

^/ / 2 N 2 TT' 2 N 2 dUAiip'^), .^s 

Sm("^m) = gM^MimM) = 9m \p^=ml, 

and where mjv/ is the meson mass. At this point we take the mesons to be spinless for the sake of simplicity. The 
generalization to mesons (or baryons) with arbitrary spin is straightforward. 




FIG. 1: Diagram describing the meson mass operator. 
To clarify the physical meaning of the compositeness condition, we first want to remind the reader that the renor- 

1 /2 . 

malization constant Zj^j can also be interpreted as the matrix element between the physical and the corresponding 
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bare state. For Zm = it then follows that the physical state does not contain the bare one and is therefore de- 
scribed as a bound state. The interaction Lagrangian Eq. ([T]) and the corresponding free Lagrangian describe both 
the constituents (quarks) and the physical particles (hadrons) which are bound states of the constituents. As a result 
of the interaction, the physical particle is dressed, i.e. its mass and wave function have to be renormalized. The 
condition Zm = also effectively excludes the constituent degrees of freedom from the space of physical states and 
thereby guarantees that there will be no double counting. The constituents exist in virtual states only. One of the 
corollaries of the compositeness condition is the absence of a direct interaction of the dressed charged particle with the 
electromagnetic field. Taking into account both the tree-level diagram and the diagrams with self-energy insertions 
into the external legs (that is the tree-level diagram times Zm — 1) yields a common factor Zm which is set to zero. 
This allows for another interpretation of the compositeness condition in as much as the condition Zm = leads to the 
correct normalization of the electric form factor of a charged particle at zero momentum transfer. By using the Ward 
identity which relates the electromagnetic vertex function at zero momentum transfer to the derivative of the mass 
operator for the on-mass-shell particle and taking into account the compositeness condition in the form of Eq. ([5|), 
one obtains 

A-(p,p) ^ ^^^f^ = ^v^''^^ = 2p^. (6) 

A^(p,p) is the zero momentum transfer electromagnetic vertex function. The form ([6]) is more suitable for analytical 
calculations because the derivative of the mass operator with respect to determines the electromagnetic vertex 
function at zero momentum transfer. 

B. Inclusion of photons 

The interaction with the electromagnetic field is introduced in two stages. The free Lagrangian of quarks and 
hadrons is gauged in the standard manner by using minimal substitution: 

9^Af± ^ (9^ T ieA'')M^, d'^q (9^ - ze,A'')g, a^g -> (9^ + ieqA^'^ (7) 

where e is the positron (or proton) charge and where is the quark's charge (e„ = | e, = — ^ e, etc.). Minimal 
substitution gives us the first piece of the electromagnetic interaction Lagrangian 

Cf'^W = Y.^,A^{x)J^{x)+eA^{x)J';^{x)+^A\x)M-{x)M+{x), 

J^ix) = q{x)-f''q{x), Jlj{x)=i(Ad-{x)d''M+{x)-M+{x)d^'M-{x)y (8) 

It is important to reiterate that there is no direct coupling of the photon to the meson in our relativistic quark model 
due to the compositeness condition Zm = as has been emphasized in Sec. HIl A. 

Gauging the nonlocal piece of the Lagrangian in Eq. ([T]) proceeds in a way suggested in P^] . In order to guarantee 
local gauge invariance of the strong interaction Lagrangian, one multiplies each quark field q[xi) in £f*j with a gauge 
field exponential according to 

^str+em(2)^^^ = ffMM(x) j dx^ JdX2FM {x , X,, X2)qi{xi) g-.i ^(^^ .-^^) Fm Am e"*^"^ ^(^^■^■^^ 52(^1), (9) 

where 

I{x„x,P) = J dz^A^(z). (10) 

X 

It is readily seen that the full Lagrangian is invariant under the local gauge transformations 

q,{x) ^ e*'='''^("''g,(a;), g,(x) ^ g,(a;)e~*^'''^("'', M{x) ^ e"'"'^'^''^ M{x) , 

Af'ix) ^ A^'{x)+^^'f{x), (11) 

where cm — &q2 ~ 6gi is the meson electric charge. 
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The second term of the electromagnetic interaction Lagrangian arises when one expands the gauge exponential 
in powers of up to the order of perturbation theory that one is considering. Superficially the results appear to 
depend on the path P which connects the end-points in the path integral in Eq (fTO|) . However, one needs to know 
only derivatives of the path integrals when doing the perturbative expansion. One can make use of the formalism 
developed in J^] which is based on the path-independent definition of the derivative of I{x, y, P): 

dx^—I[x,y,P) = lim [Iix + dx,y,P')-Iix,y,P)] (12) 

where the path P' is obtained from P by shifting the end-point x by dx. Use of the definition ([T2l) leads to the key 
rule 



d 
dxf^ 



Iix,y,P)^A^ix) (13) 



which states that the derivative of the path integral I{x, y, P) does not depend on the path P originally used in 
the definition. The non-minimal substitution Q is therefore completely equivalent to the minimal prescription as is 
evident from the identities ([T2l) or ([T3| . The method of deriving Feynman rules for the non-local coupling of hadrons 
to photons and quarks was already developed in Refs. [sl and will be discussed in the next section where we apply 
the formalism to the physical processes considered in this paper. 

For example, expanding the Lagrangian Eq. Q up to the first order in one obtains (/ — wipi + W2P2) 

£™*-^^(a;) = gMM{x) dxi dx2 j dy Elj{x,xi,X2,y) A^{y) qi{xi)T M)^Mq2{x2) , 



E';,{x,x„X2,y) = / ^1^1 _i^e^Pi(^.-)-P.(^.-)+^.(.-)i5A.(p,,p,,,), (14) 



E'^ipi,P2,q) = -eg.wiiwiqf' + 21^) dW^ {-t{wiq + if - {I - t)P) 



1 

2 /I ^\l2\ 





1 



+ e,,W2(u'2g^-2r) / dt(^';,i{-t{w2q~lf -{l-t)l^) , 



Let us emphasize that the vector Ward-Takahashi identity for an off-shell photon is also satisfied in this approach: 

<Z^A^(p,p') = EM(p')-SM(y'), (15) 

where A'^(p,p') is the meson electromagnetic vertex function described by the diagrams shown in Fig. [2l In order to 
guarantee the universality of neutral and charged current transitions (conserved vector current hypothesis) one also 
has to gauge the Lagrangian with regard to the electroweak interactions. This has been described in detail in Ref. Q 
(see, also Sec. IIIip. 



C. Infrared confinement 



Let us reiterate that the relativistic constituent quark model as described up to this point has a limited range of 
applications due to the threshold constraint mi -f 7712 > rriM- In particular, processes involving light vector mesons, 
e.g. the p and K*, and excited states cannot be treated within the model. In this section we extend the applicability 
of the relativistic constituent quark model by taking into account quark confinement effects. In order to elucidate how 
confinement is implemented in our framework we begin by considering a scalar one-pole propagator. The Schwinger 
representation of the propagator reads 

00 

— 5 ^= / da cxpf— Q;(m^ — p^)] (16) 

TO — P J 



where the integration over the Schwinger parameter runs from to (X3. Instead of integrating to infinity we introduce an 
upper integration limit 1/A^. We call the dimensional parameter A (with mass dimension [m]) the infrared confinement 
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FIG. 2: Diagrams describing meson electromagnetic vertex function. 



scale. By means of the cutoff one obtains an entire function which can be interpreted as a confined propagator, 



daeM-a{m--p')] = ' " ^'^^^^ ^'^^''^ • (17) 



Similar ideas have also been pursued in Refs. |11| where an infrared cutoff had been introduced in the context of a 
Nambu-Jona-Lasinio model. Note that the propagator for a particle in a constant self-dual field has a form similar to 
Eq. (|17p . Such vacuum gluon configurations have been studied in and were then employed e.g. in [l3| to construct 
a model with confined constituent quarks. The propagator in Eq. (|17|) does not have any singularities in the finite 
-plane, thus indicating the absence of a single quark in the asymptotic space of states. 

However, the use of confined propagators in the form of entire functions has its own difficulties. The convolution 
of entire functions leads to a rapid growth of physical matrix elements once the hadron masses and energies of the 
reaction have been fixed. The numerical results become very sensitive to changes of the input parameters which 
requires extreme fine-tuning. For these reasons, we suggest to proceed in the following way. 

Let us consider a general /-loop Feynman diagram with n propagators. One writes, using again the Schwinger 
parameterization, 

°° n 

n(pi,...,p„) = j d^a y[d4fc]'<I>exph^a,(TO2_p2)] (^8) 

where $ stands for the numerator product of propagators and vertex functions. After doing the loop integrations one 
obtains 

oo 

n = yd"aF(ai,...,a„), (19) 



where F stands for the whole structure of a given diagram. The set of Schwinger parameters can be turned into a 
simplex by introducing an additional i-integration via the identity 

1=1 dt5{t-^a,) (20) 


leading to 

n= /d"a(5(l-^a,) F(toi,...,to„). (21) 
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As in Eq. p7)) we cut off the upper integration at and obtain 

1/A= 1 



n^= / dte-' I d"aS(l-^a,] F{tai,...,tan) (22) 







By introducing the infrared cutoff one has removed all possible thresholds in the quark loop diagram. We take the 
cutoff parameter A to be the same in all physical processes. 

In order to make contact with recent ideas on the holographic description of particle interactions we change the 
integration variable t in Eq. (|2ip to z with t = z'^. One can then interpret z as an extra space coordinate and the 
upper integration limit zir, = 1/A as the infrared scale where quarks are confined and hadronized. The analogy to 
the extra dimension holographic coordinate z introduced in holographic model is now apparent. We mention that the 
holographic model of hadrons is motivated by the anti-de Sitter/conformal field theory (AdS/CFT) correspondence, 
see Refs. There the truncation over the holographic coordinate z is necessary in order to break conformal 

invariance and to incorporate confinement in the infrared region. 

As a further illustration of the infrared confinement effect relevant to the applications in this paper we consider the 
case of a scalar one-loop two-point function. One has 



"'^^"^ J ^2 [rn^ + + ^PEr][m^ + {kE ~ ^Pe)^] ^^^^ 

where the numerator factor e"^'^^ comes from the product of nonlocal vertex form factors of Gaussian form; kE, Pe 
are Euclidean momenta. Doing the loop integration one obtains 

oo 1 

n2(p') = J dt-^j^ J da exp{ - t[m^ - a(l - a)/] + ^ - ^) V} • (24) 



The integral 112 (p^) can be seen to have a branch point at p^ = 4to^. By introducing a cutoff in the t-integration one 
obtains 

l/A^ 1 

n2(p')= / '^ \,lt)2 j ^«exp{-^[m^-a(l-a)p^] + ^(a-^)V}, (25) 



where the one-loop two-point function IV^^P^) longer has a branch point at p^ — 4m^. 

Such a confinement scenario can be realized with only minor changes in our approach by shifting the upper t- 
integration limit from infinity to 1/A^. The confinement scenario also allows to include all possible resonance states 
in our calculations. First calculations done in this paper show that the limited set of adjustable parameters of the 
model (size parameters, constituent quark masses and the confinement scale A) leads to a consistent description of a 
large number of low energy mesonic processes. We envisage a multitude of further applications as e.g. in the baryon 
sector. 



III. BASIC PROPERTIES OF TT AND p MESONS 

In this section we discuss applications of the relativistic constituent quark model including infrared confinement to 
the decays of the tt and p mesons. We start by fitting the model parameters. We calculate the leptonic constants /tt, 
gp-y and the electromagnetic couplings gT^-y-y and gpTrj- The relevant quark model diagrams are shown in Figs. [3] and 
31 The corresponding Feynman one-loop integrals read: 
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FIG. 3: Diagrams describing (upper panel) and V — >■ 7 (lower panel) transitions. 





(a) (b) 
FIG. 4: Diagrams describing P —J- 77 and V — !■ P7 transitions. 



Up" 



Nr 



(2t 



+ /"da$U-^a)(2A:+ip)^tr[5(fc)]} 



Za = a{k + ^ p)^ + (1 — a)k'^ 



V2 J (2^) 



;,(-z„)(2fc+ip)^tr[7''5(fc)]} 



<?ll?2 



(26) 



T$^(-fc2)tr i755(fc + ip)7'^^(fc- ip)7''5(fc + ip-gi) 



3 J (27r)4i 

i-i^S{k + i p)7^5(fc - i p)7"^(fc + 5 P - 91) 



(27) 



(28) 



X tr 



(29) 



where iVc = 3 is the number of colors, = 7^(^—7^) is the left-chiral weak coupling matrix and the local propagators 
S{k) etc. are defined in Eq. (|4]). 

The meson-quark coupling constants are determined from the compositeness condition Eq. (|6|). This requires the 
evaluation of the derivative of the mass operator. For the pseudoscalar and vector mesons treated in this paper the 
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derivatives of the mass operators read 



1 „ d [ d'^k 

2p^^ dp^ "J {2n)H 



$|,(-fc2)tr 



Nc f d^k , , 
2p2 J 27r 4i ' 



7^S'i(fc + wip) wi^Si{k + wip)^^S2{k ~ W2p) 



+ [nil o TO2), 



(30) 



nV(/) = 



2p2^ 7 (27r)4i 



7^5*1 (fc + wip)Y S2[k — W2p) 



— ^ /t — r:r*v(-fc )tr 



7^5i(fc + wip) wi^Si{k + wip)^^ S2{k — W2p) 



+ (mi O 777,2). 



(31) 



Because of the the p^-derivative Eqs. (1501) and (PT|) contain three propagator factors. The evaluation of the one-loop 
integrals Eqs. (P5)) -(I5T |) proceeds as described in Appendix A. 

There are four adjustable parameters: the constituent quark mass 771, = mu(d), the size parameters K-^ and Ap, 
and the scale parameter A characterizing the infrared confinement. A least square fit to the observables yields the fit 
parameters 



A, 



A, 



A 



0.217 0.711 0.295 0.181 GcV 
In Table U we compare the results of the fit with available experimental data. 



(32) 



TABLE I: Basic properties of the tt and p meson. 



Quantity 


Our 


Data [151 


/,., MeV 


130.2 


130.4 ±0.04 ±0.2 




0.23 


0.276 


9pi 


0.2 


0.2 


QpTT-y, GeV^""" 


0.75 


0.723 ± 0.037 



As a first application of our approach we calculate the pion electromagnetic form factor Ft^ generated by the 
diagrams in Fig.[5J and the pion transition form factor F^r-y^* generated by the diagram in Fig. UJa). In the first case, 
we are interested in the space-like region = —Q^- In the second case one photon is on- mass-shell qf = and the 
second photon has a space-like momentum squared q| = —Q^. The electromagnetic radii are related to the slope of 
form factors at the origin = — 6F'(0). Our result for the electromagnetic radius r^r = 0.612 fm is in good agreement 
with the present world average data of r^r = (0.672 ± 0.008) fm [l^. The result for the radius of the transition form 
factor r^^ = 0.315 fm^ confirms the monopole-type approximation of the CLEO data [l^ and is close to the CELLO 
measurement 17] of r^^ — 0.42 ± 0.04 fm^. 

The electromagnetic pion transition form factor is displayed in Fig. [5] and compared with data from DESY [l8|. 
the Jefferson Lab CoUaboration 19] and the CERN NA7 Collaboration In Fig. [S] we display the F^^^.{Q'^) 
form factor in the space- like region up to 4 GeV^. 

Note that the slopes of the theoretical curves are quite sensitive to variations of the size parameter Ajr . In order to 
exhibit the sensitivity to A^ we plot three curves for A^ — 0.711, 1 and 1.3 GeV in Fig. [5l Similarly, we display the 
sensitivity of the results on the infrared confinement scale for values of A = 0.010, 0.181 and 0.200 GeV. 

One should mention that our FTr^^-{Q'^) form factor behaves as 1/Q^ at large in accordance with perturbative 
QCD. The calculated form factor behavior disagrees with the new data above 4 GeV^ presented by the BABAR 
Collaboration f2ll]. For a recent theoretical analysis of the BABAR data, see e.g. Refs. p^ . 
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FIG. 5: Form factor tt^ — >■ 7r*7* as function of the space-li ke p hoton momentum . Data are taken from the JLAB _FV 
Collaboration [13], DESY ^ and CERN NA7 Collaboration [2o[. 



CLEO 
CELLO 
our result 








) . 4 r 







35 - 






) . 3 - 


> 












o 





25 - 


a 




).2 - 













15 - 








a 










) . 1 - 







05 - 



CLEO I — 
CELLO I — 
our result 




1.5 2 2.5 3 3.5 4 



FIG. 6: Form factor 77* — >■ tt" as function of the space- like photon momentum Q^. Data are taken from CLEO [l€ 



IV. EXTENSION TO OTHER MESONS AND SOME APPLICATIONS 



A. An extension to strange, charm and bottom flavors 



In this subsection we extend our approach to mesons containing strange, charm and bottom quarks. We accordingly 
have to introduce a new set of fit parameters, namely, the values of constituent quark masses rtiq (q = s, c, b) and 
the values of size parameters Am for the corresponding mesons. Note that we keep the value of the confinement 
scale parameter A fixed for all physical processes. For the fitting procedure we choose the weak and electromagnetic 
leptonic decay constants, see. Tables ITlllIVI We use the electromagnetic leptonic decay constants for the neutral vector 
mesons defined by fv — rnygv-y where the dimensionless constant gv-f was introduced in Sec. Illll bv Eq. (|27p . Some 
comments should be done with respect to our fitting procedure. The values of strange, charm and bottom quark 
masses and of the size parameters of hadrons are determined by a least square fit to available experimental data for 
weak and electromagnetic leptonic decay constants shown in Tables [TT] and IIIII In addition, we include the decays 
Vc — ^ 77 and J/ip — ijcj (see Table ITVT) for a better determination of the charm parameters. Since the size parameters 
Aj/^ and A^^ are found to be close to each other we assume that the unknown value of A^^ is approximately the 
same as Ax • The description of the rj and 77' mesons is more involved because both contain a non-strange and strange 
quark-antiquark admixture which is characterized by the mixing angle 9p. Due to SU(3)-breaking we distinguish 
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between the size parameters of the non-strange and strange pieces of the quark current, i.e. 

7] — > - ^ sin ^ {uu + dd) - cos 5 $a^, ss , 
T] — > + cos (5 $A^/ (wM + dd) — sin 5 $A^/ ss , 

(5 = Op -01, 6l/ = arctan^. (33) 

We use the experimental data of the electromagnetic decays involving rj and 77' (see Table HV)) to fit the size parameters 
and the mixing angle. It appears that the best value for the mixing angle 6'p is —18°. This result is consistent with 
the range of values between —10° and —20° [l^ obtained in direct extractions of 9p from decay data involving 77 and 
77'. A mixing angle of or very close to —18° is also found in a direct analysis of the ratio r(77' — ^ 27)/r(77 — > 27) [l^ 
or of tensor meson decay widths [26j]. The least square fit yields the fit parameters: 



nib 



0.360 1.6 4.8 GeV 



^ri Arf' A^, Ax ^K* ^D* ^D^, ^D* ^B/B* ^BJB; Aj/^ \a/B^ 



(34) 



(35) 



0.295 0.70 0.85 0.27 0.45 0.87 0.30 1.4 2.3 1.95 2.6 3.35 4.4 3.3 3.0 5.07 5.0 GeV 
In Tables HIWVl we list the results of our numerical fit to the weak and electromagnetic leptonic decay constants 



TABLE II: Weak leptonic decay constants fp{v) in MeV. 



Meson 


Our 


Data fl5] 


7r~ 


130.2 


130.4 ± 0.04 ± 0.2 


K' 


155.4 


155.5 ±0.2 ±0.8 


D+ 


206.2 


205.8 ±8.9 


Dt 


273.7 


273 ± 10 


B~ 


216.4 


216 ±22 


B'l 


250.2 


253 ± 8 ± 7 


Be 


485.2 


489 ± 5 ± 3 [23] 


P+ 


209.3 


210.5 ±0.6 fl5] 


D*+ 


187.0 


245 ± 20l 2 [24] 


DI+ 


213 


272 ± 161^0 [^ 


3*+ 


210.6 


196 ± 24+ f [24] 


Bf 


264.6 


229 ± 20l 16 [24] 



together with their experimental values. Once the fit parameters are fixed one can use them to calculate a wide range 
of electromagnetic and dilepton decay widths. The results of the calculation are presented in Table IIVI which also 
includes experimental results whenever they are available. 

Note that the electromagnetic properties of the tt, K, p and uj mesons (form factors and radiative transitions) have 
also been considered in the light-front constituent quark model (LFCQM) developed in Ref. [l^,!!!]. It was shown 
that a reasonable description of data can be achieved with the use of the following values of the constituent quark 
masses niu — rud ~ 220 MeV and nis = 410 MeV, which are close to our values ?ti„ — nid = 217 MeV and = 360 
MeV. 



B. Dalitz decays 

In this section we apply our approach to the Dalitz decays P — )■ and V — Pl~^l~ (for a theoretical review, 

see e.g. (29|). In particular, we analyze the transition form factors of the Dalitz decays rj — > -fii'^ii~ and uj — > n'^ fj,~^fi~ 
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TABLE III: Electromagnetic leptonic decay constants fv of vector mesons with hidden flavor in MeV. 



Meson 


Our 


Data [15] 




148.0 


154.7 ± 0.7 




51.7 


45.8 ± 0.8 




76.3 


76 ± 1.2 




277.4 


277.6 ± 4 


T(ls) 


238.4 


238.5 ± 5.5 



decays. Both form factors have been measured the SERPUKHOV-134 Collaboration (Protvino) [1^ and recently by 
the NA60 experiment at the CERN SPS [30]. The ui — > 7r"/i+/i^ transition form factor has also been analyzed by the 
SND CoUaboration at the BINP (Novosibirk) 

The diagrams describing the Dalitz decays are shown in Fig. [71 They include both the diagrams with direct 
emission of the photon and resonance diagrams with an intermediate V ^ ^ transition. Note, that such decays have 
been studied in {32j by using a similar relativistic quark model with confinement (Quark Confinement Model). The 




FIG. 7: Diagrams describing the Dalitz decays P -fl^l (upper panel) and V — > Pl^l (lower panel). 

differential cross sections w.r.t. the dilepton mass squared — {pi+ + Pi-Y reads (the data is usually plotted w.r.t. 
the dilepton mass M = ^-S- 0) 



2ar(P-^77) 



3 TT 



Fp{q^ 



(' + 5^)(-^)"('-^ 



(36) 



dr{v Pi+i- 



laT{V 
3^ 



P7) 



/ 2m? \ / 4mr\ ^/^ / 



1/2, 



{niv — mp)' 



3/2 



3/2 



(rny + mpY 



Fviq') 



4nif < q^ < {my — nip)'^ 



(37) 



In our approach the normalized transition form factors are calculated in terms of the diagrams Fig. [T] The form 
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TABLE IV: Electromagnetic and leptonic decay widths in keV. 



Process 


Our 


Data [15] 


tt" — > 77 


5.40 X 10"^ 


(7.7 ±0.6) X 10^^ 


V ^ 11 


0.51 


0.510 ±0.026 


77' 77 


4.27 


4.28 ±0.19 


Vc 11 


4.55 


7.2 ±0.7 ±2.0 


Vb -> 77 


0.43 




p° ^ e+e- 


6.33 


7.04 ± 0.06 


Lj — >■ e^e~ 


0.76 


0.60 ±0.02 


(j> — !> e^e^ 


1.27 


1.27 ±0.04 


J/^p — >• e+e" 


5.54 


5.55 ±0.14 ±0.02 


T ^ e+e" 


1.34 


1.34 ±0.018 


— > 7r*7 


72.42 


68 ±7 


P° ^ ^7 


63.25 


62 ± 17 


W — >■ 71^7 


682.71 


703 ± 25 


oj — >■ r/7 


7.63 


6.1 ± 2.5 


77' —5- oj'y 


12.44 


9.06 ± 2.87 


-> t;7 


51.72 


58.9 ± 0.5 ± 2.4 


</) — > ?7'7 


0.41 


0.27 ±0.01 




40.86 


50 ±5 




122.04 


116 ± 10 




0.62 


1.5 ±0.8 




20.27 


< 0.9 X 10^ 


d:± ^ i:»?7 


0.30 


< 1.8 X 10^ 


—.4- 

5*=^ ^ 5*7 


0.36 




B*° B°y 


0.12 




Bf ^ 


0.12 




J/ip ^ rid 


1.89 


1.58 ± 0.37 


T ^> T]bl 


0.02 





factors are given by 

Fp{q^) = X [gp^^iq^] + V gpv^{q^) Dv{q^) gv^{q^)] , (38) 



V=p, cu, ^ 

Fv{q^) ^ ^ y< {gvPjiq^) + 9VPV'{q^)Dv'{q^)q^gv'yiq^)}- (39) 

^^^^ V'=p,LO,<t> 

The decay constants gp^^ = gp-y-yiO) etc., are the same as defined in Sec. IIIII One has to note that if wc assume 
that the quantities g{q^) in Eqs. ((38|) and ((39)) do not depend on q^ then we reproduce the resuhs of vector meson 
dominance (VMD). 

The contributions of the intermediate vector meson resonances are described by a Breit-Wigner form where we 
have used a constant width in the imaginary part, i.e. we have used 

Dv{q') = — j^- ^ . (40) 

In Fig. [S] we plot the transition form factors for the two Dahtz decays r] — )■ /i+/x^7 and lu — > Tr^/i+zi^ measured by the 
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SERPUKHOV-134 Collaboration (Protvino) [H, the CERN NA60 Collaboration ^ and the SND Cohaboration 
BINP (Novosibirk) j31|]. For completeness we also display our prediction for the Dalitz decay rj' 76+6" in com- 
parison to the older data given in [2^. In Table IVl we present our results for the slope parameters defined by F^{0). 




FIG. 8: The calculated form factors of the Dalitz decays ui tt^I'^1~ , rj 'yl'^l~ and r/' — > ■yl^l' as functions of the dilepton 
mass M = Experimental data are taken from [23. [30t for uj — s- n^'n'^fi^ and 77 — !> 7/1+^", from [sj] for r/ — >■ 7e"'"e~ and 

from [29!] for rj' — >■ 7e"'"e~ (the five uncorrected by background points are included). For comparison we plot the VMD-curves. 

Finally, our predictions for the Dalitz decay widths are given in Table IVll The agreement with the data is generally 
quite good except for the decay tt" — )■ je'^e~ which comes out too small in our model. This is not unexpected since 
the decay 7r° — >■ 77 also comes out too small in our model (see Table IV). 



TABLE V: Slope parameters of the Dalitz transition form factors in GeV 



Decay mode 


Our 


Data 




1.4 


5.5 ±1.6 [33] 


rj 7Z+P 


1.4 


3 ± 1 [34] 
1.95 ±0.17 ±0.05 [30] 


V ^ ji+r 


1.1 


1.68 [35] 




1.9 






1.9 


2.24 ± 0.06 ± 0.02 [22,] 


^ T]i+r 


2.2 






2.2 






1.4 






2.5 





V. SUMMARY 



We have refined a relativistic constituent quark model developed in our previous papers to include quark confinement 
effects. Quark confinement was implemented by introducing an upper cutoff on a scale integration which, in the original 
quark model, extends to infinity. The introduction of such an infrared cutoff removes all physical quark thresholds 
in the original quark diagrams. The cutoff parameter is taken to be the same for all physical processes. We adjust 
the model parameters by fitting the calculated quantities of the basic physical processes to available experimental 
data. As an application, we calculate the electromagnetic form factors of pion and the transition form factors of the 
Dalitz decays P — and V Pl^l^ . We extend our approach to mesons containing strange, charm and bottom 
quarks and calculate their leptonic and radiative decay constants. 
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TABLE VI: Dalitz decay widths in keV. 



Decay mode 


Our 






Data [15] 


7r° 7e^e~ 


6.4 X 


10" 


-5 


(9.39 ±0.72) X lO""' 


rj — )■ 7e+e~ 


8.5 X 


10" 


-3 


(8.84 ±1.14) X 10"^ 


V -> 7M^M" 


0.4 X 


10" 


-3 


(0.40 ±0.06) X 10"^ 


rj' — )■ 7e^e^ 


9.0 X 


10" 


-2 


< 0.2 




2.0 X 


10" 


-2 


(3.2 ±1.2) X 10"^ 




0.66 






< 2.3 




7.0 X 


10" 


-2 




a; — >■ TT^e+e^ 


6.23 






6.54 ±0.77 


+ - 


0.65 






0.82 ±0.20 


p° — >■ Tye^e" 


0.46 






< 1 


uj rje^ e~ 


5.6 X 


10" 


-2 


< 9.3 X 10-2 


4> iqe^ e~ 


0.44 






0.49 ± 0.04 


(j) ri'e^e^ 


2.1 X 


10" 


-3 






2.2 X 


10" 


-2 


< 4 X 10^2 
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Appendix A: Loop integration techniques 

In order to demonstrate how the loop integrations are done consider a n-point one-loop diagram with n local 
propagators Si{k + Vi) and n Gaussian vertex functions $i ± Wj+n)'^)- In Minkowski space the one-loop diagram 
can be written as 

--2-tr J|$, (-(fc±t;,+„)2) TiSi[k + Vi) (Al) 

i—l 

where the vectors Vi are linear combinations of the external momenta pi to be specified in the following, k is the loop 
momentum, and the F^ are Dirac matrices (or strings of Dirac matrices) for the ith meson. The external momenta 

n 

Pi are all chosen to be ingoing such that one has ^ = 0. Due to translational invariance the integral Eq. (jAip is 

1=1 

invariant under a shift of the loop momentum fc — >■ fc + Z by any four- vector /. The four- vector / may be any linear 
combination of the external momenta pi . 

Using the Schwinger representation of the local quark propagator one has 

00 

S^{k + Vi) = {m,+ jt+j^,)x Jd^, exp[-P,{m^ ~ (k + v^f)] . (A2) 



For the vertex functions one takes the Gaussian form. One has 

(-(fc±w,+„)2) = exp [/3j+„ (fc + u,+„)2] i = l,...,n, (A3) 
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where the parameters f3i+n = Si = 1/Af are related to the size parameters. The numerator factors rrii +j/i+'i!fi can be 
replaced by a differential operator in the following manner: 

/ 74 7 " f ? " 1 

— tr J] / d(3, e-^-'"'(m,+ ^+ x exp <^ ^ ft (fc + ^;,)' 

* i=li [i=l J 



/flit. " °r / i\ f 2" 

-2- tr n / rfA e-^'"' f m,+ + - ^ J x exp <^ /3fc2 + 2fcr + ^ ft 

i— 1 Q 2—1 

2n 2n 

„2 , 



(A4) 



where P = ^ Pi and r = Z^* '^i- Next one does the loop integration and moves the Gaussian exp I — r^//3 j to the 

1=1 i=l ^ ^ 

left of the differential operator. By using the identity 

-^ + E/3''^'=« E ft (A5) 

one obtains 

oo ( 

^ f d8 ' I ^ 1 

/„(pi, ...,p„) = n / ^ ^ ~ E ft""*' + « E ft ft - ^j)' 

i=ln I j=l l<j<j< 



<j<2 n 



tr[|r, (m,+ ) . (A6) 



X 

i=l 



We have written a FORM jSg program that achieves the necessary commutations in a very efficient way. As described 
in the main text the set of Schwinger parameters ft (i — 1, ...,n) can be turned into a simplex by writing 

oo oo / " \ 

y"d"/3F(ft,...,/3„) = y"dir-iy"rf"a<5( 1-^a, J F(t a„). (A7) 

\ 1=1 / 

One then arrives at the final representation of the one-loop n-point diagram in the form 



/„(pi, ...,p„) ^ dt 



i=l 

nrJm,+ ^,--l^ ^,,V (A8) 



n 

X tr 

i=l 



where 



2ioc = E caajAij, 

i—l l<i<j<n 



Zl 



Z2 = 



n 2n 

E E ft ~ E ' 

i—l l<i<j<n 



ft ft- , 

n+l<i<j<2 n 



n 2n 

r = t ^ tti + S ^ ft 

■i—l i— n+1 
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Here, 



Si/S, 



i=l 



The matrix A, 



''^j)^ (1 ^ hi l£ 2r7,) depends on the invariant variables of the 



process. 

There are ahogether n numerical integrations, n — 1 a-parameter integrations and the integration over the scale 
parameter t. For the derivative of the (two-point function) mass operator one has to do one more a-parameter 
integration due to the extra propagator which comes in after the differentiation. The integration of the derivative of 
the mass operator proceeds in analogy to the n point function case described in this Appendix. We mention that the 
correctness of the numerical integration procedure can be checked very conveniently by shifting the loop integration 
momentum by a fixed momentum four-vector. The numerical evaluations have been done by a numerical program 
written in the FORTRAN code. 

Some further remarks are in order. The convergence of the loop integral Eq. (|A8p is defined by the local a form 
zioc- If zioc < the ^-integration becomes divergent due to contributions from the large t-region. The large t-region 
corresponds to that region where the singularities of the diagram with its local quark propagators appear. However, 
as described before, if one introduces an infrared cutoff on the upper limit of the t-integration, all singularities 
vanish because the integral is now convergent for any value of the set of kinematical variables. Note that the one-loop 
integration techniques described in this appendix can be extended to an arbitrary number of loops in a straightforward 
manner. Of particular interest is the extension to the two-loop case needed for the description of baryon transitions. 



Appendix B: Gauge invariance of the p° f transition 



In this appendix we want to demonstrate by an explicit calculation that the transition amplitude p*^ — >■ 7 written 
down in Eq. (1271) satisfies gauge invariance, i.e. that e.g. the on-shell photon has only two transverse degrees of 
freedom. According to the formalism developed in Sec. |TT]B there are the two contributions to the transition — > 7. 
The first contribution results from the minimal substitution in the free quark Lagrangian and is depicted in Fig. Ela) . 
In Fig. [2Ib) we depict the point interaction contribution resulting from gauging the nonlocal Lagrangian. 

We begin by considering the transition of an on-shell p to an off-shell photon with invariant mass p^. The corre- 
sponding transition amplitude M^'^{p) must satisfy the gauge invariance condition p,yM'^''(p) = as has already been 
assumed in writing down Eq. (j27|) . We shall now show by explicit calculation that the non-gauge invariant pieces 
in the two contributions cancel each other resulting in an overall gauge invariant contribution. First we isolate the 
non-gauge invariant pieces in the two respective contributions by writing 



(£k_ 



<I>,(-fc2)tr(7^^(fc + ip)7''5(fc-ip)) 
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The non-gauge invariant contributions are given by 

00 

dt 



i^'Kp') 



(Bl) 



e , zi = tm 



st p^ 



dtt 



{s + tf 




dae 



1 



Z2 — t ^m^ — a{\ — a)p^^ 
whereas the gauge invariant contribution is given by 



st f 1 



s + t\ 2 



(B2) 



Itip') 



dtt 







dae 



1 2t^ f I 
a — — 



2 (s + i)2 V 2 



(B3) 
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Similarly we have 



oo 1 

dt s 



1 



- 1- 



2 sa\ sa 



s+t V s+tj s+t 4 



da ( 



1/ 2sq;\sq! 
1 - 



4 V s + t J s + t 



Z3 = im — ( 1 — 



s a \ s a 



s + tJ 4 



(B4) 



Here s — l/A^. The non-gauge invariant pieces li^''^^ and 7^'^'' cancel each other as can be seen by the following 

transformations. First, we note that the integrands of the integrals li^\p'^) and l'f'\p'^) may be expressed via the 
derivatives of Z2 and Z3, respectively. 



I'^\p' 



2t^ 



s + t {s + tf 



S + t 



\ - [a 



1 



\dz2 
2) da 



li'Hp') 



1 P^ f ^ 2as\ as 



s + t 4V s + t J s + t s + t 



1 - a 



dz3 
da 



The a-integration can be done by using the bound-state conditions Z2{a = 1) = Z2{a = 0) = Z3(a = 1) = zi. One 
obtains 



00 1 

dtt 



{s + 1)3 



dtt 



(s + t)3 



(3)/„2n 



00 1 

dt s 



(s + t)3 



a e 



dt 



{s + t)3 



Finally, one has 



ii'^{p') + e\p') + ij,'Hp')^ 



dt 



(s+tr 



1 - 



s+t s+t 



0. 



(B5) 



The gauge invariance condition pi,M^''(p) — also guarantees that the longitudinal component of the photon decouples 
as — > 0. 
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